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Abstract. Wc study certain generalized Cauchy integral formulas for gradients 
of solutions to second order divergence form elliptic systems, which appeared in 
recent work by P. Auscher and A. Rosen. These are constructed through func- 
tional calculus and are in general beyond the scope of singular integrals. More pre- 
cisely, we establish such Cauchy formulas for solutions u with gradient in weighted 
L2(R-4^"', t"dtdx) also in the case |a| < 1. In the end point cases a = ±1, we show 
how to apply Carleson duality results by T. Hytonen and A. Rosen to establish 
such Cauchy formulas. 



1. Introduction 

A fundamental problem in modern harmonic analysis has been if the Cauchy sin- 
gular integral on a Lipschitz curve defines an L2 bounded operator. Calderon [8] 
showed that this is indeed the case when the Lipschitz constant is small, and Coif- 
man, Mcintosh and Meyer [9] showed boundedness for any Lipschitz curve. In 
connection with the latter work it was also realized that boundedness of the Cauchy 
integral, a problem in harmonic analysis, was equivalent to the Kato square root 
problem, a problem in operator theory, in one dimension. The relation can be 
seen as follows. The Cauchy singular integral on the graph of a Lipschitz function 
(7 : R — 7- R is given by 

u{y){l + i9'{y))dy 



p.v.- 

TT 



sgn{B D)u{x) 



u e L2(R). 



/R {y + w{y)) - {x + ig{x)) 

The operator theoretic expression sgn{BD) for the Cauchy integral on the right hand 
side is interpreted as follows. In L2(R) we have the self-adjoint differential operator 



D :-- 



' dx-' 



or equivalently the Fourier multiplier — and the accretive multiplication 



operator B = {l + ig'{x))~^. This yields a bisectorial operator BD which was shown 
to have a bounded holomorphic functional calculus, see Mcintosh and Qian [T3] . 
In particular the bounded symbol sgn(A) = ±1, ±ReA > 0, yields an L2-bounded 
Cauchy integral operator sgn{BD). Note that when B = I, this formula is simply 
the Fourier relation ;^J-'(p.v.l/x) = sgn(^). 
The Kato square root estimate 



d 



d 



■—a[x)—u 
dx dx 





du 


2 


dx 



on the other hand follows from the boundedness of sgn{BD) for more general accre- 
tive coefficients B. In higher dimension the Kato square root estimate || divAVM||2 

||Vm||2 follows from a similar estimate ||sgn(i?D)|| < 00, with B = ^. and 



^ Formerly Andreas Axelsson. 
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D 



div 
-V 



A major difficulty in higher dimension, n > 2, is that D has an 

d 



infinite dimensional null space. Note that 



^ dx 

dx 



i-^ when n = 1, if we iden- 



tify ranges = C. In higher dimension, the Kato square root problem on R" was 
solved by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian [H] and the more 
general result that operators of the form BD have bounded holomorphic functional 
calculi was proved by Axelsson, Keith and Mcintosh [7j. 

Coming back to the Cauchy integral, in this paper we study certain generalized 
Cauchy type operators which appeared in recent work by Auscher and Axelsson [1] . 
More precisely, the aim is on one hand to give some complementary results for these 
Cauchy operators on certain weighted L2-space between the end point cases studied 
in [T], and on the other hand to show duality results in these end point cases, using 
results of Hytonen and Rosen [TT] . 

Our Cauchy operators are constructed in the above spirit, by applying suitable 
bounded and holomorphic symbols to an underlaying differential operator like BD. 
We shall even need to apply more general operator valued symbols to the differential 
operator which, changing the setup slightly, will be of the form DB = B~^{BD)B. 

To formulate the problem, consider a divergence form second order elliptic system 

divt,xM't^ x)^t,xU = 

in the upper half space R^"^*^ := ; t > 0,x G R"}, n > 1. We as- 

sume that u : R^"^"" — >■ C", m > 1, is vector valued and that coefficients A G 
Loo(R+*'"; £(C'^^"'""''"*)) are accretive in the sense that there exists k > such that 

inf Re{A{t,x)v,v)/\v\'^ > K 

t)GC(i+")'"\{0} 

for almost every (t, x) e R++'^. With minor modifications, all our results are valid 
under a weaker Carding type inequality, uniformly in t. See [1]. 
A natural gradient of solutions u is the conormal gradient 



ViiM 



where di,^u = {A'Vt,xu)± denotes the conormal derivative and VyU = VxU denotes 
the tangential gradient of u. Similarly, divy = diVa; and curly = cur^ will denote 
tangential divergence and curl. We write and v^^ for the parts of a vector v normal 
and tangential to the boundary. 

Question 1.1. For solutions to a given divergence form equation divt^xAit, x)Vt,xU = 
as above, is there a Cauchy type formula 

for the conormal gradient? 

To answer this question, we first need to specify function spaces for VaU- We shall 
use the following natural subspaces of L]^'^(R}^^). Here and below, we often suppress 
the range of functions in notation, for example L2(R+'*'") = L2(R+'''"; C^^"*""^™). We 
also write || ■ ||l2(R") = II ' lb- 
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Definition 1.2. Define 



{/ : R^+" ^ C(i+")™ ; \\NiW,f)h < oo}, 



using the non-tangential maximal function Nf{x) := sup\y_^\^g\f{s,y)\ and L2 

Wliitney averages W2f{t, x) := t~^^~^'^^/'^\\f\\L2(^wit,x)) over Wliitney regions W{t, x) : = 
{{s,y) ; 1/2 < s/t < 2, \y-x\<t}. 



For -!<«<!, let 

L2[ 



{/ : R^+" ^ C(^+")"^ yy^^ < 00}. 



It was shown in [H Lem. 5.3] that 



(1) 



1 

sup - 

<l+n\ 



\\fs\\lds<\\NiW2f)\\ 



2 < 



,ds 



s\\2 



We think of A''2,2(R+ ) as a substitute for L2(R+ , t"") in the endpoint case a = —1, 
which allows for non-zero traces. 

To state our results, we next introduce the operators that we use. For more details, 
see [Ij. With the second order divergence form operator divt,xA{t,x)'Vt,x comes a 



first order self-adjoint differential operator D :- 







diV|| 





parallel to the boundary R", and a pointwise transformed coefficient matrix 



acting tangentially. 



(2) 



B 



a 

ca~ 



1 d 



ca 



if A 



a b 
c d 



How these operators appear are explained in Section [2l They act on C'-^"'""^'"- valued 
functions, written as column vectors with normal parts first and tangential parts 
second. We write ft{x) = fit,x) for such functions in R^"^*^, and similarly for the 
coefficients Bt{x) = B(t,x). Write 

£t{x) = £{t, x):=I- Boixr'B{t, x), 

where Bo{x) are some t-independent accretive coefficients Bq E Loo(R"; /^(C^^^"^™")). 
Often Bo{x) = B{0,x), but not always. 

Our fundamental operator is DBq. Both as an operator in L2(R") and in L2(R++", t"), 
acting in the x- variable for each fixed t > 0, it defines a closed and densely defined 
operator with spectrum contained in a bisector S^^ = 5'^^+ U (—5'^+), where 

5-^+ := {A e C ; |argA| < U {0}, w < 7r/2. 

In [7j it was proved that DBq has a bounded holomorphic functional calculus, which 
gives estimates of operators h{DBQ) formed by applying holomorphic functions h : 
S*^ — )■ C, a; < /i, to the operator DBq. In particular, we shall need the operators 



A : 



\DB^\ 



e-*^,t > 0, 



X±{DBo), 

Pt POO 

Sft := / Ae-^'-'^^'E^fsds + / Ae^^'^'^'^E^ f^ds. 
Jo Jt 



For the first three operators, we view DBq as an operator in L2(R") and apply 



the scalar holomorphic functions A H- |A| 



±A, ±ReA > 0, A ^ e"*!^! and 



A I— )■ x±('^) 1 if i Re A > and elsewhere. For the definition of S, we view DBq 
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as an operator in L2(R+'^'^, t") and apply the operator- valued holomorphic function 
A I—)- F{\), where 



When |a| < 1, F is a bounded function, and hence S = F{DBq) is a bounded 
operator on L2(R+^",t°). For the endpoint spaces a = ±1, see Section [21 For 
further details of this operational calculus, see [U Sec. 6]. 

Theorem 1.3. Consider first the case \a\ < 1 and assume that divt,xA{t, x)'Vt,xU = 
with estimates ||/||i2(R^+" f^) < ^ of the conormal gradient f := VaW- Then there 
exists a function G Eq L2(Ii.^) , such that 

(3) ft = A'^e-*^Eo+/i+ + SStft, a := (a + l)/2. 

From this follows estimates 



l'2t 

(4) supt-M \\h~''fs\\lds< 

t>0 Jt 



L2(R.i7",i") 



and we have limits 

|2, 



(5) lim ||A-"/, -/i||^cis = 0= limrW ||A^"/,||^ds 

■ 0+ i^oo 



where h := h}~"e~'^^EQ£sfsds, \\h\\2 < II/IIl2(R++" If furthermore a > 0, 

then we have the pointwise L2(R") estimates sup^^o l|A^'^/t||2 ^ II/IIl2(r^+" t") '^'^^ 
limits 

(6) hm ||A-/, - /.||2 = = hm ||A-V*||2. 

Conversely, ^/ ||^^||l^(r,i+") sufficiently small, then the Cauchy type formula 

(7) := (/ - S£r'A'^e-*''E+h+, h+ e E+L^iRn, 

constructs a function f with \\f\\L2{R^^"t") ~ ll^^lh? which is the conormal gradient 
f = of a solution u to (\iVt^x^{tiX)'W t,xU = 0. 

When a = +1, the above holds with the following changes. We need to assume 
throughout that \\S\\^ < oo, and for the converse statement that \\S\\^ is sufficiently 
small, where || ■ ||* denotes the Carleson-Dahlberg norm from Definition \2.1[ Here 
0" = 1 and we have traces in L2(R") sense. 

When a = —1, the above holds with the following changes. We need to assume 
throughout that < oo, and for the converse statement that is sufficiently 
small. Furthermore we need to replace L2(Ri'(_'^"', t") throughout by A''2,2(R+'")- Here 
(T = and we have traces only in the square Dini sense ^ . 

Note that the trace spaces for / are exactly the fractional homogeneous Sobolev 
spaces iJ~°'(R"). We record the following result, proved in Section [31 

Proposition 1.4. Let < a < 1. For all f E D(A~°") with curln/n = 0, we have 

l|A~VI|2^||/b-.. 
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Note that in the cases a = ±1, the t-independent coefficients Bq are uniquely 
determined by B, see [U Lem. 2.2]. When |a| < 1, this is not the case. 

We also remark that the representation formula ([3]) can be used to prove various 
other estimates of solutions. See [1]. 

Example 1.5. To recognize (JTj) as a Cauchy formula, consider the special case 
A = B = I,n = l= m and a = —1. Then Vu will be anti-analytic, and hence 



is the Cauchy type reproducing formula we look for. Letting £ = 0, a = and 
B = I in Theorem II. 3[ formula (I7j) reduces to 



e~*'^'x-i-(~'^); ^iid that D is the Fourier multiplier — ^. 

As less trivial example, we consider the special case m = 1, A being real t- 
independent coefficients and a = —1. Then it was shown in [H] that 



is the conormal gradient of the single layer potential, for normal vector / scalar fields 
h e L2(R"). Here r(s y) denotes the fundamental solution to div t,xA{x)'Vt,x in R^^" 
with pole at {s,y). Note that in the case of the Laplace equation A = I, 



is the Cauchy/Riesz kernel, (T„ denoting the area of the unit sphere in R^"*"". 

Finally, we remark that for general systems, m > 2, and general coefficients 
A, the operators defined from DBq by functional calculus are usually beyond the 
scope of singular integrals. For example, the known constructions and estimates 
of the fundamental solution T(^s.y) require De Giorgi-Nash local Holder estimates of 
solutions to the divergence form equation, which may fail for systems, m > 2. 

The end point cases a = ±1 and the estimate of S in the case |a| < 1 was proved 
in pp. In this paper we supply the details of the remaining results stated for |a| < 1 
in Section [H] and a simplified proof of the estimate for S in the case a = +1 in 
Section [21 In the final Section HI we make some remarks on applications to the 
Neumann and Dirichlet problem for divergence form equations. 



The aim with this section is to give a simplified proof of the estimates of the 
singular integral operator S = Sa from [T] in the case a = +1, using Carleson 
duality results from [TT]. We start by deriving the integral equation Q for the 
conormal gradient / = Vau from the divergence form second order differential 
equation diYt^xA{t, x)Vt,xU = for the potential u. 

Splitting A as in ([2]) , we have /x = adtU+bV i\U and /y = VyW. Thus the divergence 
form equation, in terms of / reads 






1 

(T„, (t2 + |a;|2)(H-n)/2 



2. Carleson estimates of operators 



dth + div,{ca-\U-bf,) + df,) = 0. 
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The condition that / is the conormal gradient of a function u, determined up to 
constants, we express as the curl-free condition 

'9,/||=V||(a-i(/x-fe/ii)), 
curly /ii = 0. 



In vector notation, we have 

dt 



7x" 




diV||" 






+ 


-V|| 





-1 



a —a 
ca^^ d — ca^^b 





7x" 




./ii. 



0, 



together with the constraint curly /y = 0, or in short hand notation 

dtft + DBtft = 0, fte R{D) =: n. 
With the t- independent coefficients Bq, we rewrite the equation as 
(8) djt + DBoft = DBoSJt. 

We shall use freely known properties of operators of the form DBq. See [21 I]. 
In particular DBq is a (non-injective if n > 2) bisectorial operator and L2(R") = 
HiDBo)®^,, where 

We now integrate the vector- valued ordinary differential equation ([H]) for ft & T-L. 
Applying the projections Eq, we have 

dJt--Afr = -AE^SJu 
where ft^ ■= E^ ft. 

Formally, assuming lim(_^o+ ft = /o and lim^^oo ft = 0, we integrate these two 
equations 

and subtraction yields the integral equation 



-iA 



E^f, + ^^i/i. 



In Section Owe show by a rigorous argument that, depending on the function space 
for /, integration indeed yields this equation with Eq f^ = K^h^ , for some G 
£"(^172 (R"). In this section, we discuss estimates of the singular integral operator S 
and the multiplier in particular in the case a = +1. 



Definition 2.1. For functions / in R^"^", define the Carleson functional 



Cf{x) := supr " / / \g{s,y)\dsdy 

r>0 J J\y-x\<r-s 

\f{s,y)\s~'^dsdy, x G R". Define the 



y—x\<s 



and the area functional Af(x) 
Banach space 

fRi+") := {/ : R^^'^ ^ 0^^+")'" ; ||C(W^2/)||l,(R" 



2.21 



P(W^2/)|U.(R")<00}. 
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The equivalence of the Carleson and area functionals 

||C^||l^(R") ~ \\Ag\\L,in"), K p < oo, 

follows from [lOl Thm. 3]. From pT] Thm. 3.2, 3.1], we recall the following duality 
result. 

Proposition 2.2. The Banach space A^2,2(R'+'") is the dual space o/ 6*2,2 (R'+''*) ™- 
der the L2(R+"^") pairing. The space C2 2(R-+'^") is not reflexive, that is N2 2(R+^")* ^ 

By Proposition 14.51 and ([I]) we have continuous inclusions 

L2(R^+",t-^) C iV2,2(R++") and C2,2(R++") C L2(R^+",t). 

The estimates for the singular integral 5* are as follows. 

Theorem 2.3. // |a| < 1, then S : L2(R++",r) ^ L2(R++",t") is a bounded 
operator. For a = ±1, we have bounded operators S : L2(R+'^", t""*^) — ?■ A''2,2(R+'") 
and S : 6*2,2 (R++") ^ i^2(R++", t). 

This result was proved in [1], with the exception that the Carleson space C2,2(R+^") 
was not known there and in the endpoint case a = +1 only the estimate 

was proved. We here survey the proof from [l] and supply the missing estimates of 
S and S separately in the case a = 1. 

Proof. First recall the rigorous definition of the singular integral 5* from 
[H Sec. 6,7]. For fixed e > 0, define truncated singular integral operators 

/•£ roc 

SJt:= / vt{t,s)Ae-^'-'^^E^fJs+ / v7 (t, s)Ae~^-'-'^^E^ f^ds, 
Jo Jt 

where r]f are compactly supported approximations of the characteristic functions 
of the triangles {{t,s) ; < s < t} and {{t,s) ; < t < s}. Then Se : 
L\°'^(R_,_; L2(R")) -> L^(R+; L2(R")) is a well defined operator. More precisely, 
we define ?7°(t) to be the piecewise linear continuous function with support [l,C)o), 
which equals 1 on (2, 00) and is linear on (1, 2). Then let ?7e(t) := ?7°(t/e)(l — ?7°(2et)) 
and rif{t, s) := ri^{±{t - s) /e)r]^{t)r]^{s). 

For |a| < 1, it was proved in P, Thm. 6.5] that 5*^ are uniformly bounded and 
converge strongly to an operator S in L2{^^^, t"). The idea of proof was to view 5*^ 
as being constructed from the underlaying operator DBq in L2(R+'^'^, t") by applying 
the operator- valued symbol A H- -Fe(A), where 



Fe{\)ft:= / ry+(t,s)|A|e-(*-^)l^lx+(A)/,ds+ / r^:{t,s)\\\e~^^~'^\>^\x-{\)fsds, 
Jo Jt 

yielding = F^{DBq). It was shown by Schur estimates that 

sup ||i^e(A)||i2(R^+",tc)^L2(R3^+",ta) < 00, 

for any u < fi < 7^/2. Moreover, for any fixed e > there is decay 
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and F,{\)f -> F{\)f for each A G and / G L2(R++",t") as e 0. As shown 
in [H Sec. 6.1], from this and square function estimates for DBq it follows that 
= F^{DBq) are uniformly bounded operators in L2(R+^",t") which converge 
strongly to ^ = F{DBo). 

At the end point space a = —1, the above bounds fail on L2(R}^^,t~^) and we 
split the operator as 

S, = F^{DBo) + FXDBo), 

where 

F'Wft := [\t{t,s)\\\e~^'-^^\\^i\)fjs 
Jo 

/oo 

rt+2e 

- / iv,{t)v,{s) - v:{t,s))\\\e-^^^'^\\.{\)fsds 
Jo 

and 

POO 

F^{X)f, := r/,(t)e-*l^l / Veis)\X\e-'\\-Wfsds. 

Jo 

On one hand, the term F^ can be treated as in the case \a\ < 1, yielding a bounded 
operator F^{DBq) : L2(R++",t"^) L2(R++", t"^)- On the other hand, the term 
F factorizes bounded operator 

L2(R^+",ri) ^ L2(R") ^ Ar2,2(R^+"), 

where the second factor does not converge strongly in £(L2(R"), iV2,2(R++")) but 
only in £(L2(R"'), L2{a, b; L2(R"')) for fixed but arbitrary < a < 6 < oo. 

Now finally consider the end point space a = +1. Here the appropriate splitting 

is 

S, = F!{DBo) + F^{DBo), 

where 

F!Wft := rr/+(t,s)|A|(e-(*-)l^l -e"(*+^)l^l)x+(A)/.ci. 

^0 

POO 

- / {rie{tMs)-r^tit,s))\X\e-'^'+'^\\AX)fJs 

Jt-2e 

/oo 
rj:{t,s)\X\e-^^-'^\\.{X)fsds 

and 

POO 

F^iX)fr.= vMMe^'^'^ / Veis)e'-^^\+iX)fsds. 

Jo 

(Note the duality F^ = {F})* and F^ = (F^)*.) On one hand, the term F^ can be 
treated as in the case |a| < 1, yielding a bounded operator F^{DBo) : L2(R+'''", t) — )■ 
L2(R+'^",t). On the other hand, the term F^ factorizes as a bounded operator 

C2,2(R^+") -> L^iRn ^ L2(R^+",t). 
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The bounds of the second factor follow directly from square function estimates, 
whereas the bounds of the first factor follow by duality from the non-tangential 
maximal estimates used in the case a = —1. We have 

(POO \ POO 

Note that in this case, both factors converge strongly as e — )■ 0, using that compactly 
supported functions in R^"^" are dense in (72,2(1^4^") • Proof of this density result is 
in [TH Lem. 2.5]. Thus S, ^ S strongly as operators C2,2(Ry"") ^ i^2(R+^", t). □ 

For a multiplier S, clearly 

sup ll^/llL,(Ri+",t") = 11^ 



L2(R!+",t") - ll'-llLoo(Ri+") 



for any a. For a = ±1 we have the following more refined Carleson multiplier 
estimates. Define the Carleson-Dahlberg norm 

||£:||.:= ||Clyoo(?)||^^ 

using Loo- Whitney averages Woog{t,x) := \\g\\Loo{W{t,x))- 

Theorem 2.4. The following are equivalent. 

(i) S : f{t,x) S{t,x)f{t,x) IS bounded A/'2,2(R++") ^ L2(R++",ri). 

(ii) £ : f {t,x) £{t,x)f {t,x) IS bounded L2(R++",t) C2,2(R+"^"). 

(iii) S has the Carleson-Dahlberg estimate < 00. 

// this hold, then 

\\£\U^ sup \\£f\\L2{R\+\t-^)- sup II^/IIc2 2(rV"")- 

Proof. The equivalence of (i) and (iii) follows from [TTl Thm. 3.1]. The equivalence 
of (i) and (ii) follows from Proposition US] since S* : iV2,2(R+^") ^ L2(R++",t-^) is 
the adjoint of £ : i:2(R++", t) C2,2(R++'') under the L2(R++") pairing. □ 



3. Proof of the Cauchy formula 

In this section, we prove Theorem 11.31 The end point cases a = ±1 were proved 
in [1]. Thus it remains to show the results for |a| < 1. We remark though that 
in [T], in the case a = 1 it was only shown that S£ : L2(R+'^'^,t) — )■ L2(R+'^'^, t). 
With the intermediate Carleson space 6*2,2 (R+^") available now from [H], we have 
the refined mapping result 

£ : L2(R^+",t) ^ C2,2(R++") and S : C2,2(R++") ^ L^iKl+^t), 

if ll^^ll* < 00, proved in Section [2l 

Proof of the representation formula ([3]). Assume that |a| < 1 and that divt,xA{t, x)'Vt,xU = 
with estimates ||/||i2(R,i+" t") < 00 of the conormal gradient / := Vau. As in P, 
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proof of Thm. 8.2], for any e > 0, integration of ([H]) gives 

+ e-W {e-'^-e-^'-'^^)E+fsds-2e e-^'-'^^E^ f^ds =: I - II + III - IV, 

with equality in L2(a, 6; %) for any fixed < a < 6 < oo, where / converges to / in 
L2{a,b;7i). Using estimates ||e~*^ — e~(*~^^'^|| < s for III and ||e~'^''~*)^|| < 1 for IV, 
we obtain 

\\IIIh<e^^-''^/^^0 and ||/F||2 < e^^+^^Z' ^ 0, 

as e — )■ 0, for each fixed a < t < b. Since S^Sf converges in L2(R]^^,t°') by [H 
Thm. 6.5], it follows from the equation that II converges in ^2(0, b; %) for any fixed 
< a < 6 < 00. 



Write ft := lim^^oe ^'^e ^ J^^ Eq fgds. Since f = f — S£f by the equation, 
L2(Ri+",t-) ^ ll/llL2{RV+",i«) < Using the identity / = A JI^{sAe-'^yds/s on 
T-L, and e"*^/, = /t+s so that ||/t+s||2 < ||/s||2, we estimate 

((A*)->, /J = 4 / ((.A*)2-'^e-^^>, s'^e-^Vj- < 
io s 

for any (p G □((A*)"'"). Thus = A^'ht with supt>o ll^itlh < ||/||L2(R^+",t^)- We now 
want to let s — !■ in the identity 

Integrating against a test function G 1^2 (o, ^; "H), we get 

\^tJt+s)dt= (^j\*e-'''*Mt,h, 

By continuity of translations in ^2(0, b; "H), the left hand side converges. Since func- 
tions of the form A* e~^^* (ptdt are dense in L2(R") and hs are uniformly bounded, 
hs h weakly in L2(R") when s -> 0. We conclude that ft - SSft = ft = A'^e'^^h 

■ /'■Rl+n . 



in LorRi+",r). □ 



Proof of the estimate (jl]) . By equation ^ it suffices to estimate the weakly singular 
integral operator 

/t 1*00 
= / Ai-'^e-l^-*IX'"^*"'V.^^s + / Ai-'^e-(*-^)^(/ - e-2^^)Eo+/,ds 

Jt/2<s<2t Jo 
poo p 

+ / Ai-'^e-(^"*)^(J - e~^'^)E^fJs + e"*^ / 

J2t Jb 



R\[t/2,2i] 

=:I + II + III + IK 
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We first estimate IV by duality. For cf) G L2(R"'), we have 

\{IVA)\< / ||(sA)^"-e-^^>|b||sV.||2-<||</>||2||/IL,(Ri..^,.), 
Jo s + 

so that supi^o ll-^^lb ^ II/IIl2(R++"'*")' estimate II, we note that 

||^l-<Xg-(t-s)A(j_g-2sA)|| 

This gives 



SUp||//||2< / s/e-'^Wfshds 
t>0 Jo 



t/2 \ 







L2(R^+",t") II La (R+ 



A similar estimate applies to III. We are left with the local weakly singular integral 
I, which we estimate 



l/IU < 



Jfshds 

t/2<s<2t 



\t- S 



l-o- 



2* ds ^ 



It/2 

if a > 0. If a < 0, we at least obtain the weaker estimate 

'u/2<s<2u 



M — S| 



u/2<s<2u 1^ -^l"^ ^ J \Ju/2<s<2u 1^ -^l"*" 
Jt \Ju/2<s<2u P ~ -^l 



•4t 

't/2 \Js/2<u<2s 1^ ~ -^l 

This proves the estimate of US'/til 2 and therefore of (j4]). □ 



Jt/2 \Js/2<u<2s I" ~ -^1 / 



Proof of Theorem \1.3[ We supply the remaining arguments in the case \a\ < 1. 
Having established the representation formula ([3]), write this as 



f, = A-ie-'^E+ + S£tft). 

The estimates above for S give the stated estimates for ||A^°"/t||2- For the traces, it 
remains to prove that 

\im{S£f)t= / A'-''e-'''E^£sfsds=:h-, MmiSS f)t = 0, 
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either in L2(R") or square Dini sense. By the estabhshed uniform bounds, we may 
assume that ft^O only if a < t < 6 for some 0<a<6<oo. In this case, 

J a<s<min{t,b) J max{t,a)<s<b 

f W-^'e-^'-'^'^E+Esfshds < [\t - sy-'Uhds ^ 

J a J a 



Since 

rb 



when t — 7- oo and 

rb 



J a J a 

as t — )■ 0, we have proved the traces. 

The converse result is obtained by reversing the argument leading to the repre- 
sentation formula ([2]), outlined in Section [21 Note that square function estimates 
give 



r u^e-'^E^hmedt = r mre^'^E^hm'^ < wh^ 

Jo Jo ^ 



|2 

|25 



and if ||^^|lL2(R^+",t«)^L2(R^+",tc.) < ||£'|Il^(r^+") < 1, then I - SS is invertible on 
L2(Ri+",r). □ 



Proof of Proposition L4_. Consider the operator A* = \BqD\. From square function 
estimates and accretivity of Bq it follows that 

\\\B*D\g\\2 ^ \\B*Dgh ^ \\Dgh ~ \\\D\gh, 

where 

' ' [ (-V„div„)V2j 

Thus D(A*) = D{\D\), and interpolation shows that D((A*)'") = D(|D|'") with equiv- 
alence of norms. See remark following [4], Thm. 4.2], where interpolation between 
homogeneous norms gives the estimate || (A*)'^5f||2 ~ || |Z^|°"5'||2, for all g G D((A*)'^) = 

D(|Dr). 

By duality, we have for / G D(A~'^) n H that 
l|A^VI|2= sup {A~'^f,{A*r9) 

96D((A*)-),||(A*)-g||2 = l 

~ sup (/,^) = ||/||^_,, 

36D(|Z)|<'),|||D|'^g||2=l 

since |||D|'^(7||2 = Hs'II^/.t for (7 G "H and \D\'^g = for (7 G "H"*". □ 

4. Applications to the Neumann and Dirichlet problem 

It is important to note that in the previous sections, we have always worked with 
the quantity 



VnU 



the conormal gradient of a solution whole. On the contrary, for the Neumann 

and Dirichlet problems we need to work with the two components d,^^u, the Neumann 
datum, and VyW, the Dirichlet datum, separately. In doing so, we leave the functional 
calculus of the bisectorial operator DBq and go beyond the Cauchy integral. 
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T af-rt l+n^ 
±1 I , 



+ ^- 1 iVo.2fRi+"), a = -l. 



Consider the function spaces 

From Theorem 11.31 and Proposition 11.41 it follows that there is a well defined and 
bounded trace map 

VyiM(t,x) H- Vam(0,x) 
taking solutions u of daYt^x^{tiX)V t,xU = with Va^* G -^fl-R-i^"); ^'^ ^a^Ir^^ E 
H~'^{'R,'^). This is true for any bounded accretive coefficients A when |a| < 1, and 
in the endpoint cases a = ±1 we need to impose the Carleson-Dahlberg condition 
x)— ^4(0, < oo. If furthermore we assume smallness of a;)— yl(0, x)||oo 
in the case |a| < 1, or smallness of ||y4(t,x) — A(0,a;)||^, in the case a = ±1, then we 
have equivalence of norms 

In this case, the Hardy type subspace 

E+if"'^(R") := {VauM c ij-'^(R") 

is a closed subspace of i/~°"(R"). We make the following definitions. 

Definition 4.1. Let a G [—1, 1] and a = (a + l)/2 G [0, 1]. Consider bounded and 
accretive coefficients A(t, x), and assume smallness of A{t, x) — A{0, x) in the above 
sense depending on a. 

Let WP{Neu, H'") denote the set of coefficients A for which the Neumann map 
Ej;ii"-'^(R"; 0(1+")"^) ^ ij-'^(R"; C™) : Vam|r" ^ (VA^lRn)^ is an isomorphism, 
so that in particular it has lower bounds 

+ ~ II Va^Ir" ||H-<T(I^n) < l|C^I^A'"|R"llij-'^(R")- 

Let WP{Dir, H^~'^) denote the set of coefficients A for which the Dirichlet map 
E+H-''{R''; C(i+")™) ^ H~''(W; C^™) : Vam|r" ^ (Vam|r")|| is an isomorphism, 
so that in particular it has lower bounds 

ll^^^llLf(R^+") ~ l|VA^^|R"||i/-<T(Rn) < ||V||n|Rn||^_,(j^„). 

When A = A{x) are t-independent, the Hardy subspace E^H^'^(TV^) is the range 
of projection = x+{DBq), which acts boundedly in if~'^(R"). We have that 

L^3Bo^ E^ G /:(ij-'^(R")) 

is locally Lipschitz continuous. More generally for t-dependent coefficients, we see 
from Theorem 11.31 that 



/■oo 

Eig = E^g+ / Ae-^X"^.((J - 
Jo 



ds 



is a projection onto the Hardy space E\H '^(R") along the null space Eq H '^(R" 
In this case we have that 

_i_ _L,, ^ W^Wt i'-r1+"^i < 1, 



'A -^0 ll_H'-<^(R")-i>H-<'(R") 



a = ±1. 



Since in this way, the Hardy space of solutions depends continuously on the coeffi- 
cients, we obtain the following perturbation result. 
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Proposition 4.2. If we have a well posed boundary value problem for t -independent 
coefficients A = A{x) G WP{Neu, H~'^) , then there exists e > such that A G 
WP{Neu, H~'^) whenever sup ^-^^-^1+11 \ A{t,x) — A{x)\ < e when \a\ < 1, and when- 
ever sup^gR^n |y4(0,x) — A{x)\ < e and \\A(t,x) — A{0,x)\\^ < e when a = ±1. 
The corresponding result also holds for the Dirichlet problem. 

Example 4.3. The optimal case is when a = 0, in which case all coefficients belong 
to WP{Neu,H-^/^) and WP{Dir, H^^^). This is a simple consequence of Gauss' 
theorem, which yields 

/ / {A{t, x)Vt,xU, Vt,xu)Ve{t)dtdx 

JJr\+'- 

= 2e / / {d^^ut,ut)dxdt - e^^ / / ( 

Jl/(2e) JU" Je 7R" 



{du^ut, ut)dxdt, 

l/(2e) ■ 

with ?7e(t) as in the proof of Theorem 12 .Si Taking limits e — )■ 0, by accretivity of A 
this gives the estimate 



\\ft\\idt<\mjH-^M\ifo) 



1/2 ~ 



of the conormal gradient /. Since max(|| (/o)x||ij-i/2 1| (/o)|| ||h-i/2) ^ UoWh 
II/IIl2(R+"''")' absorb either factor of the right hand side, on the left, and 

obtain the claimed lower bounds. 

Example 4.4. Much more subtle are the endpoint cases a = ±1. The Neumann 
problem with data in L2(R'^) {a = —1) is usually denoted (A^)2 in the literature, 
and the Dirichlet problem with data in if^(R") (a = —1) is usually denoted {R)2 
and referred to as the regularity problem. See Kenig [121 Sec. 1.8]. 

The Dirichlet problem with data in L2(R") (« = +1) is usually denoted (-0)2- 
In this case our Definition 14.11 differs from the standard one in that we require the 
square function estimate 

/ / \Vu\Hdtdx < / \u\'^dx 

rather than the non-tangential maximal estimate which usually defines (-0)2- See 
Kenig [121 Sec. 1.8]. It should be clear from Section [2] why we prefer to define the 
Dirichlet problem through the square function estimate here. 

Note that by maximal function estimates proved in [H Thm. 2.4], A G WP{Dir, L2) 
implies that (-0)2 holds, modulo one subtle point. There are coefficients A G 
WP{Dir, L2) where even for good Dirichlet data (p G LiilV) n H^^^{W), the solu- 
tions G L2(R++",t) and u^/^ G L2(R++") are distinct. See [21 Sec. 5] and [H]. In 
this case, there will exist some Sobolev space < a < 1/2 where A ^ WP{Dir, H'^). 
Here u^^"^ is seen to be the solution obtained from Lax-Milgram's theorem, which 
is the one for which non-tangential maximal estimates are required in the problem 
{0)2- However, the solution which is estimated in [H Thm. 2.4] is . 

There are coefficients A for which these endpoint boundary value problems are 
not well posed. For positive results, it is known that WP{Neu, L2), WP{Dir, H^) 
and WP{Dir, L2), as well as the less well known fourth end point boundary value 
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problem WP{Neu, H contain all t-independent coefficients A{x) which are Her- 
mitean, A*{x) = A{x), or of block form, A ^ 



d 



or are constant A{x) = Aq. 



There is also a duality result for these boundary value problem, which reads as 
follows. 

Proposition 4.5. For the Dirichlet problem, we have 

A G WP{Dir, H") if and only if A* G WP{Dir, if ^~'"). 

For the Neumann problem, we have 

A G WP{Neu, H~") if and only if A* G WP{Neu, H"-^). 

This can be proved as in [5l Sec. 17.2]. We remark that it is well known that the 
Dirichlet problem (-0)2 holds whenever the regularity problem {R)2 holds, whereas 
the reverse implication is not true in general. The reason that this reverse implication 
holds in Proposition 14.51 is that we require a stronger square function estimate 
rather than a non-tangential maximal estimate for the Dirichlet problem with data 
in LsfR"). 
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